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We study a class of generalized inflation models in which the inflaton is coupled 
to the Ricci scalar by a general f(ip,R) term. The scalar power spectrum, the spec- 
tral index, the running of the spectral index, the tensor mode spectrum and a new 
consistency relation of the model are calculated. We discuss in detail the issues of 
how to diagonize the coupled perturbation equations, how to deal with an entropy- 
like source, and how to determine the initial condition by quantization. By studying 
some explicit models, we find that rich phenomena such as a blue scalar power spec- 
trum, a large running of the spectral index, and a blue tensor mode spectrum can be 
obtained. 
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1 Introduction 



Inflation [lj has been very successful in solving the problems in the standard big bang 
cosmology such as the horizon and flatness problems. Fluctuations are created during 
inflation, which provide seeds of large scale structure and the Cosmic Microwave Back- 
ground (CMB) anisotropies having been observed by COBE and WMAP satellites, 
the latter will be measured more precisely by PLANCK. 

In spite of the remarkable achievements, there are still open problems, such as a 
natural realization of inflation models in a fundamental theory, and how to further 
pin down a specific model including the number of the inflatons and the form of the 
potential. In ref. [2], some theoretical problems are listed. For example, the hierarchy 
of scales, the trans-Planck problem, the singularity problem and the candidate of the 
inflaton. As for the data fitting, there is still a significant power deficit at I — 2, and if 
the running of the spectral index is introduced as a parameter, the WMAP three year 
data mildly favors a blue spectrum with a large running [3], which is not common in 
inflation models. In the past, a few models were proposed to accommodate a large 
negative running such as the noncommutative inflation model [3]. For other models 
we refer to [5], El [7J E] • Discussions on models violating the null energy condition and 
leading to a blue tensor spectrum and possibly a blue scalar spectrum, as well as on 
the WMAP three year constraints, can be found in [9] 

To have a broad scope, and to be prepared for possible surprises in the future 
experiments, it is essential to study classes of models with novel physics built in, for 
example, quantum gravity and string theory may induce higher derivative corrections. 
In this paper, our primary goal is to carry out a more careful study of a class of models 
in which in addition to the Hilbert-Einstein action and a minimal coupled inflaton, 
there is a general nominimal coupling of the form f(ip, R). Some of these models were 
studied before. However, when we examined the literature carefully, we found that 
only some special models were reliably studied. These include: f(R) gravity models 
[TT] and the f(<p)R gravity models [12] . a unified analysis of these two cases were made 
in [13] . Furthermore, even for f(R) action and f(<p)R action, no concrete models were 
investigated in detail. Remarkably, the result for the tensor modes power spectrum in 
[T3] is applicable to more general cases. We note, however, that the pioneering work 
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on inflation with an additional scalar curvature squared term was done long time ago 

m 

In [14] . the models with a nonminimal coupling term —^f(^R)if 2 are studied, and 
the calculations are carried out when / is quadratic in R and the nonminimal coupling 
gives dominate contribution in the rolling of (p. It is shown that a blue spectrum with 
a running larger than that of a common minimal model can be produced. However, 
the constraints on the fluctuation of the inflaton and the scalar fluctuations in the 
metric are over-simplified in [T^] . 

In this paper, we study a general class of models in which the coupling between 
the inflaton <p and the Ricci scalar R is described by a general function f((p, R). Both 
the background dynamics and the perturbations are studied, and the power spectrum, 
the spectral index, the running of the spectral index and the perturbations for the 
tensor modes are calculated, which can be related to the experiments today and in 
the future. 

Some ingredients in the calculations are standard as being used in the minimal 
inflation models ( [15] is a nice review). Due to the higher derivative terms, the con- 
straints among the fluctuations are no longer purely algebraic. We develop a general 
method to deal with this problem. Diagonalization of the perturbation equations 
is no longer simple. Upon assuming a slow-roll inflation, we can still have a nearly- 
conserved quantity which reduces to the usual curvature perturbation when the higher 
derivative terms are suppressed. Thus, although there is an entropy-like perturbation, 
we can bypass this problem by using this new nearly-conserved quantity. Also, since 
we can not directly write down the action for the perturbation modes, we provide 
an argument in relating the perturbation in a diagonalized perturbation equation to 
the canonically normalized mode, thus set the initial quantization condition. Some 
of these methods are quite general and can be used in other generalizations of the 
minimal inflation models. 

By calculating the power spectrum and related quantities, we find that there are 
rich phenomena in this class of models. We can have a blue spectrum and a large 
running of the spectral index in some models, while in other models the spectrum 
stay red and the running can become positive. Also, in principle, one can have an 
initial fast-rolling period when the higher derivative terms are comparable to the 
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Hilbert-Einstein term. The consistency condition between r and n-r is modified. A 
blue tensor mode power spectrum can appear even in some small inflaton models 
in some parameter region while in minimal inflation models the tensor spectrum is 
always red. 

The paper is organized as follows. In Section 2, we derive the equations of motion 
of background quantities, and spell out the conditions of slow-roll approximation. In 
Section 3, we write down the equations for perturbations. After diagonalizing and 
simplifying the equations to obtain a Bessel equation, we propose a nearly-conserved 
quantity which has the form of comoving curvature perturbation in the minimal 
inflation model limit. In Section 4, we study the initial conditions by relating the 
fluctuation under study to the canonically normalized mode. Since the canonical 
variables are normalized by quantization, this procedure help us to fix the integration 
constants in the solutions of the Bessel equation. In Section 5, the power spectrum, 
the spectral index, the running of the spectral index, and the ratio of the tensor and 
scalar perturbations are given. In Section 6, we calculate several classes of models 
and discuss the properties of these models. We conclude in Section 7. 

2 The Model and Background Dynamics 

We consider the following action with nonminimal coupling: 



We consider first the evolution of the background, a homogeneous and isotropic 
spacetime with a space-independent inflation field. We work with the Friedmann- 
Robertso n- Walker (FRW) metric 




(1) 




(2) 



where t is the comoving time and r is the conformal time. 

It follows from the action ([!]) that, the equations of motion are 



(3) 




(4) 
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where p and p are energy density and pressure respectively, as the diagonal elements 
of the spacial averaged energy-momentum tensor, and can be written as 

M 2 1 RF — f 

P = -fC^ + V + ^l-SHF), (5) 

P = ^(^ 2 -V-^ T ± + F + 2HF), (6) 

where F = -J^f(ip, R), and dot denotes the derivative with respect to t. 
We shall impose the following slow-roll conditions: 

H 

6 S -A<1, (8) 

7 - (9) 

The first and second slow-roll conditions are standard as in the minimal inflation 
model. The third slow-roll condition is new in the nonminimal model. It is a natural 
requirement, because F can be written as a function of R and ip and should be slow- 
rolling during inflation. Moreover, it will be shown that 7 appears in the spectral 
index and the ratio of tensor and scalar perturbations. To be consistent with the 
CMB data, it should be small in any viable candidate inflation model. In addition to 
the above slow-roll conditions, we further require that the above slow-roll parameters 
themselves are rolling slowly. 



3 Classical Evolution of Perturbations 

To calculate the scalar fluctuations, we need to perturb the Einstein equations. 
Through out the calculations, we will work in the longitudinal gauge, 

ds 2 = a 2 (-(1 + 2<f))dT 2 + (1 - 2$)dx i dx i ) , (10) 

and the fluctuation of the inflaton ip is denoted by Sip. 
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The perturbation equations are 

3#W + ^)-lvV = -^P, (11) 

3(H<f> + ip) = ^(p+p)v, (12) 
5F 

*P-<P = jr, (13) 

$ + 3H(H<f> + ij) + H<p + 2H<P + ^V 2 (0 - V) = ^Sp, (14) 

which follow from the components G^GpGUi ^ j), G\ of the linearized Einstein 
equations for perturbations, respectively. The dp, 5p and v are perturbative quantities 
in the energy momentum tensor defined as, 

T ° = -(p + 5p), 7? = -(p + p)v ti , Ti = 3(p + 5p), (15) 

with the explicit form 

M 2 1 1 

5p = ^[<p50+~(-f^ + 2V ip )6<p-3H5F + (3H + 3H 2 + -V 2 )5F 

+ (3#F-0 2 )0 + 3F(#0 + ?/O], (16) 

M 2 1 .... 9 

6p = -J-[ip5ip+-(f^-2V ><p )S<p + 6F + 2H6F+(-H-3H 2 -—V 2 )6F 

J7 Zi OCX 

-F4> - (<p 2 + 2F + 2HF)<P - 2F(H(j) + ip)], (17) 
M 2 

(p + p)v = -^(<p5<p + 5F - H5F - F0) (18) 

By using equations (1121) and ( ITBl to cancel Sip and terms in (flTi) respectively, 
one gets 

F(0 + ^) + (ifF + 3F - — )(0 + ^) 

+ [(3#F + 3F) - —(HF + 2F) - ^V 2 U 

<fi a z 

+[UHF + HF-F)- —(HF — F) — ^V 2 U = (19) 

ip cr 

And by using 5F = Fr5R + F jV 8<p and F = F^i? + F >v >ip, we can write ({TBI) in a 
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more explicit form 



F(2HF - HF + F)(4>-ip) 

+ [6F, R (4HH + H)- F][F(<j> + </>) + (HF + 2F)<P + (HF - F)^} 
2F >R (2HF -HF + F)x 

[3i> + 3#0 + 12Hj) + Q{2H 2 + H)<j> + \v 2 (<p - 2^)] 



(20) 



In the limit of minimal inflation model, (120!) reduces to the constraint equation ip = <f>. 

Equations f)19p and (I20p are coupled differential equations for ip and 0. In some 
special cases, for example, F = F(cp) or F = F(R), these equations can be diag- 
onalized to a single differential equation [13J. While in a more general case which 
is of main interest in this paper, we have to develop a new method to solve these 
equations. 

To diagonize the equations, we define 



where ipk and 4>k are Fourier components of ip an d 0, and the explicit fc-dependence 
of Qk will be shown later in this section. We can Taylor expand Qk as 



where Go and G\ stand for the zeroth and first coefficients in the Taylor expansion. 

So under the condition that the nonminimal coupling is not too strong com- 
pared with the minimal coupling gravity, the above expansion works well near the 
horizon-crossing and we will see in the following sections that Go will be the dominant 
contribution. It ought to be true for physically interesting models where the energy 
scale of inflation is smaller than the Planck scale or the string scale. While when the 
expansion does not work well and the high order terms of k become important, the 
high order derivative corrections will play a much more important role at the Hub- 
ble scale. In this case the dispersion relation will be modified significantly when the 
perturbation was created, and we need to modify the usually canonical quantization. 
A quantum gravity theory perhaps is in demand in order to fully understand this 
regime. 



Qk = ipk/<Pk, 



(21) 




(22) 
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We now derive Gk as a function of background quantities from (fl9l) and (j20|) . 
Substitute (I2TI) into (fl9l) and (I20I) to eliminate then use ( fl9l) and (1201) to eliminate 
the 0" terms to get a first-order differential equation. Differentiate this equation once 
and subtract with ( fl9l) . we can get another first-order differential equation. From 
the consistence condition of the two first-order differential equations, we can get the 
expression of Q k , which to the lowest order is the solution of the following equation: 



{[FF + 24HF, R {HF - HF)\g k + F(F - 24MF fi )} 

x{[2HHF(F - \2H 2 F yR ) - AHF yR (AHF - 5HF)^ - FF^}Q k 

k 2 k 2 

-2H(F - 12H 2 F R )(HF - HF) + 8HF R (HF + HF) — - FF—} 

a 2 a 2 

= {[2HF(F - 12H 2 F R ) - 2F R (2HF - HF)^]G k 

k 2 

-2(F - 12H 2 F R )(HF - HF) - 2F R (2HF - HF) — } 2 , (23) 

a 2 

we see that indeed Q k is a function of only the background. By solving Go and G\ 
from this equation, it can be shown that Go and G\ are slow-roll quantities and Go 
has the simple form 

G = 1 - ™ (24) 
FH 

With the use of the linear relation between ip k and fa, to leading order of slow-roll 
parameters, we derive a differential equation for a single variable, say <p k 

h " k + {aH) I H(i + g k ) + 26 + 37 ) ^ + e(pk 
^ I h4t^ + mtt-& + T& + »V-° (25) 



By introducing a new variable 



the equation ( 1231) can be rewritten as 



_ F*(i + g k ) A 

u k = ; fa, (26) 



1 f r , 3 — Q k AQ k \ 2 



r 2 V 2 (1 + Qk) l + Gk 



where (aH) 2 has been replaced by 1/r 2 up to the leading order precision. 
Inserting the Taylor expansion of Qk, one finds 

* + h{ i+ w^-£k c h + ~** =0 - (28) 

where k is a rescaling of k, due to the contribution of the /c-dependent part of Qk, 
P 1 , ( 4G 4(1 + 2G ) \ F rH 2 

V 2 = l + {wtw 1 " "o+g^v 1 (29) 

It can be shown in the models under consideration and to the desired precision, this 
rescaling can be neglected, and we shall use k instead of k. If one wish to calculate to 
second order in slow-roll parameters in this nonminimal model, this rescaling should 
be taken into account. This effect may also arise in models other than ([1]). 
The equation (128]) can be further simplified to 

u "k + ttt ("4e + 25 - 7) u k + k 2 u k = 0, (30) 

Which takes the form of a Bessel equation. The solution of (1301) can be written as 

u k = CiV=?HW(-hr) + C 2 V^Hi 2 \~kr), (31) 

where C\ and C 2 are the integration constants to be determined by the initial condi- 
tions, and v is 

v = I v /l-2(-4e + 2 ( 5-7)^ \ - + 25 - 7 ) (32) 

There are rich phenomenona as preheating and reheating for the perturbations 
from crossing the horizon to reentering the horizon. In order not to track the whole 
evolution process, we need to find a conserved quantity well outside the horizon. In 
the minimal inflation models, the comoving curvature perturbation is introduced to 
deal with this problem. But in the models with entropy perturbation, the comoving 
curvature perturbation and its direct generalizations may not be conserved. In this 
paper, we develop a general method to derive a nearly conserved quantity, which can 
come back to the form of comoving curvature perturbation in the limit where the 
physical wave-length is large enough to ignore nonminimal coupling corrections. 

In our model, the nearly conserved quantity takes the form 

+ + (33) 
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It can be examined that lZk/(HlZk) ~ ^(e 2 ). The detailed derivation of this 
quantity is presented in Appendix A. The issues on entropy perturbation and the dif- 
ferences between the comoving curvature perturbation and the newly-derived quantity 
are discussed in Appendix B. Since the nearly conserved quantity turns smoothly into 
the comoving curvature perturbation at late time during inflation, the observables like 
the temperature fluctuations on the large-angular scales can be directly related to this 
quantity. The amplitude of entropy perturbation automatically decays well outside 
the horizon, so it does not conflict with the experiments on the observed limits of 
entropy perturbations. 



4 Quantum Theory of Perturbations 

Quantization plays an essential role in creating the perturbations from an initial 
"vacuum" state, and sets the initial conditions for the solutions to the perturbation 
equations. To do quantization, we have to expand the action up to second order in 
the perturbation variables: 

hS=\j (V) 2 - (d lV ) 2 + ^v 2 ^ d 4 x (34) 

As usual the commutation relations on the constant conformal time hypersurface are 
satified by the canonical coordinate v and its conjugate momentum. 
The equation of motion of v takes the form: 

v" k + k 2 v k + m 2 v k = (35) 

At large k/(aH), provided that the perturbation is inside the horizon and the 
nonminimal coupling effect has not been so strong to spoil the local flat space quantum 
field theory, the normalization can be determined as in the flat spacetime as 

* = ir *' T (36) 



Now recall u in the pervious section, we shall determine the role of u in this 
quantum system. Note that u can be expressed in linear combination of v and v' as: 

u k = Av' k + Bv k (37) 
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In general, A and B are fc-dependent. But up to a time-independent overall 
factor, they should be almost ^-independent. This is because the ratio of A and B 
is determined by dynamics, and during the slow-roll inflation, each mode has almost 
the same dynamics in its production and evolution. 

Insert f l37|) into fl30l) . and use fl35l) to simplify the expression. By collecting the 
k 2 v term in the resulting equation, we get A' = 0, i.e. A is independent of time. Then 
we can determine A at very late time during inflation. By comparing the formula of 
nonminimal inflation model 

Vk H 3- — v k (38) 



F§(l + G ) " Ft{l + G 0) 

with the corresponding formula of the minimal model at k/ (aH) <C 1, A is determined 
to be 

A=% (39, 

The validity of this equation can be checked numerically. In the models discussed 
in the following sections, at a few e-folds before the end of inflation, the energy scale 
is indeed low enough where the non-minimal coupling models dynamically reduces to 
the minimal model. 

Having A in hand, we are able to fix the initial condition of u. When k/(aH) is 
large, Av' k term in ([37]) is the main contribution. So the short-wave limit of u reads 

Uk = T7k e ~ lkT ^ 

kV 2k 



By expanding the Bessel function at short-wave limit, one can fix the integration 
constants in ( 13TT) 

K*l = ^£. C^O (4!) 

5 The Power Spectrum and the Spectral Index 

We now have all the pieces needed to calculate the power spectrum of the nearly- 
conserved quantity 1Z. To the leading order in slow-roll parameters, 1Z can be ex- 
pressed as 
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H H 2 

7l k ~ T k + -rrfa (42) 
H H 



M p H 2 f-kr 



-u+l/2 



We calculate 71 at a few e-folds outside the horizon, where the Bessel function solution 
still holds and 1Z can be taken as a conserved quantity from then on. The coefficients 
in the Bessel equation are calculated at k — aH. 
The power spectrum of TZ takes the form 



It follows that the spectral index is 



n s — 1 = 



= -4e + 25-7, (44) 



dliak 

As usual, n s — 1 can be calculated in another way. We can consider the explicit k 
dependence of Vn a few e-folds after crossing the horizon and do not use the horizon 
crossing condition. Then the spectral index can be expressed by v in 



n s - 1 = 1 - 2v = -Ac + 25 - 7, (45) 

which agrees with (14*4|) . It is not only a self-consistency test for the conservation of 
1Z, but also a test for the quantization. Note that the only condition for (145]) to hold 
is that the classical evolution equation (1301) remains the same outside the horizon. 
So in parameter regions where the nonminimal coupling is very strong and we can 
not expand Q as equation (1221) inside the horizon, the power spectrum fomula may 
change up to a multiplying constant, but the spectral index (j44D and the running of 
the spectral index (I46p should still be correct. 
The running of the spectral index is 



dn s 

a. = 



4e 25 7 



dink 

For the tensor modes, the results in [13] can be directly generalized to our models. 
The spectrum of the tensor modes reads 

2 H 2 

Pt = - 2lF , (47) 
12 



where we use the normalization such that in the limit of a minimal inflation model, 
the ratio of the tensor and scalar type perturbation reduces to r = 16e. 
The spectral index of the tensor modes is 



n T = _2e - 7 (48) 

In the above equation, e is always positive, so the tensor power spectrum is always 
red in the minimal inflation models. As to be demonstrated in the next section, 7 
is negative in some cases thus it can help to produce a blue spectrum of the tensor 
mode perturbation. This is a new effect in the nonminimal models. 

The ratio of the amplitude of tensor fluctuations to scalar fluctuations is calculated 
to be 

r = ^ = 8 ^ (2£+7) (49) 

The consistency relation has been modified to be 

r = ~8^n T (50) 

p 

This new consistency relation is subject to test in future experiments. 



6 Some Nonminimal Inflation Models 

In this section, we shall examine some concrete models of the following form 

f = M 2 p R + g&)R n * (51) 

During inflation, the Hubble scale is much lower than the Planck scale, the coupling of 
the inflaton with very large power of R is unlikely to be important. We will consider 
the cases when = 1,2 and study the 113 = 2 case in more detail, the latter has not 
been studied in the literature except in [14] as far as we know, and is the most natural 
candidate to make a large correction to the minimal inflation model. Of course other 
choices are interesting in their own right, although we are not going to investigate 
them in this paper. 

For the inflation potential V(ip), we consider the following potentials for the in- 
flaton: the quadratic potential V((p) = ^m 2 (p 2 , the quartic potential V(cp) = |A<^ 4 , 
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the power-law potential V(ip) = XM^ exp(— ap/M p ), the small field quadratic po- 
tential V(<p) = AMp(l — ap 2 /Mp), and the small field quartic potentialV^) = 
XMp(l — c(p A /Mp). Again, the results can be generalized directly to other kinds 
of potentials of <p. 

In our calculations, the co-moving wave number k and the nonminimal coupling 
constant Ao are taken as either input parameters or abscissas in diagrams, while the 
coupling constant in the inflaton potential V((p) is determined by COBE normaliza- 
tion 8 2 H = ±V n «4x 1(T 10 . 

6.1 Models with / = M 2 R + X (p 2 R 

As we know, this class of models can be related to a minimal inflation model with 
a single inflaton by a conformal transformation, and this kind of coupling is not of 
particular interest in quantum gravity theories such as string theory. But to compare 
with the models listed below, we shall still discuss some features of this class of 
models. 

Figfj] and FigfS] present n s — 1, a s , and r as functions of the nonminimal 
coupling constant respectively. These models generally produce a red power spectrum, 
and a very small running. The power-law potential is an exception, which is to be 
discussed in more detail in the next subsection. 

6.2 Models with / = M 2 R + \ M p 2 (p 2 R 2 

By COBE normalization condition and choosing the solution which smoothly reduces 
to the minimal case in the weak nonminimal coupling limit, n s — 1, a s , and Ut, t are 
shown in Figj3]and FigJH respectively. A red power spectrum is obtained in this case. 
A large running can be produced for some potentials when the nonminimal coupling 
is large. For the small field potential the large running is negative while the quadratic 
potential gives a positive large running. 

We note from these figures that the small field inflation with V(<p) — A(l — ap 2 ) 
is not favored by experiment when Ao — > because n s — 1 is too small and r is too 
large. But when Ao is taken into consideration, this model behaves much better in 
both figures. We also observe that ny in the small field inflation can be positive, 
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Figure 1: n s — 1 and a s as functions of Ao in the models with / = M 2 R + \Q(p 2 R. Red 
spectrum is common in this case except that the power-law potential gives a very small 
positive n s —l. Running of spectral indices are quite small except in the small field quadratic 
inflation. 
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Figure 2: and r as functions of Aq in the models with / = M 2 R + Xo(f 2 R. 
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Figure 3: n s — 1 and a s as functions of Ao in the models with / = M 2 R + \qM~ 2 (p 2 R 2 for 
the solution which has a smooth minimal coupling limit. The power spectrum is also red 
in all cases we consider. Large running can be produced when the absolute value of Ao is 
large, but only the small field quadratic potential produces a negative large running. 



0.06 ■ 

0.04 - 

0.02 - 

0.00 - 

-0.02 - 

-0.04 - 

-0.06 - 

-0.08 - 

-0.10 - 

-0.12 - 

-0.14 - 

-0.16 - 

-0.18 - 

-0.20 - 

-0.22 - 

-0.24 ■ 

-0.26 t 




— Quadratic 
- Quartic 

Power-law 

— Small Field Quadratic 



-100000 -80000 -60000 -40000 -20000 

The Non-minimal Coupling Constants 



0.75 
0.70 
0.65 
0.60 
0.55 
0.50 
0.45 
0.40 
0.35 
0.30 
0.25 
0.20 
0.15 
0.10 
0.05 



Quadratic 




Quartic 




Power-law 




Small Field Quadratic 





-100000 -80000 -60000 -40000 -20000 

The Non-minimal Coupling Constant \ 



Figure 4: tit and r as functions of Ao in the models with / = M^R + XoM~ 2 ip 2 R 2 for the 
solution which has a smooth minimal coupling limit. The small field quadratic potential 
model provide a positive ut- Such phenomena never appears in the minimal models. 
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Figure 5: n s — 1 and a s as functions of 1n(k/k c ) in the models with / = M 2 R+ XqM~ 2 (p 2 R 2 
and the quadratic potential V(ip) = ^m 2 (p 2 for the solution which has a smooth minimal 
coupling limit. Ao is taken -50000 in both figures. k c is where the COBE normalization is 
taken. We assume at k c , there are 60 e- folds of inflation. The power spectrum is red, with 
a positive running. 

while in the minimal inflation model it is always negative. If in future experiments a 
blue spectrum of gravitational waves is detected, it may be a signal that nonminimal 
coupling effect should be considered. 

Figl5] presents n s — 1 and function of ln(A;/A; c ), where k c is where COBE 

normalization is taken, we assume it is 60 e-folds from the end of inflation. To know 
the exact e-folding number at COBE normalization point requires the knowledge of 
the details of reheating, which will not be considered in this paper. 

We find that something new happens here, if we take Ao ~ —60000 and draw a 
figure like Figl5j F crosses zero at the largest scale that can be observed today, and 
the perturbations blow up during the zero-crossing. This region can not be treated 
exactly by the approach in this paper because the slow-roll conditions (Q and the 
expansion ( l22l) break down. The exact equations of motion for the background shows 
that this region is indeed unusual. This phenomenon can also happen in some other 
models studied in this paper. This phenomenon may have something to do with the 
lower CMB anisotropy power on the largest angular scales. To make sure of this, 
a detailed calculation is required. Similar phenomenon in the quintessential models 
of dark energy is studied [16] . but as far as we know, no consideration in inflation 



17 



- Quadratic 
Quartic 



-30000 -25000 -20000 -15000 -10000 

The Non-minimal Coupling Constant A, 



- Quadratic 
Quartic 



-30000 -25000 -20000 -1 5000 -1 0000 -5000 

The Non-minimal Coupling Constant k Q 



Figure 6: n s — 1 and a s as functions of Ao in the models with / = M 2 R + XoMp 2 ip 2 R 2 
for the new solution in the nonminimal model. A blue power spectrum is produced in this 
case, with a negative running when the absolute value of Ao is not large. 

models and effects on the power spectrum is made. 

The COBE normalization equation does not always have an unique solution. In 
some models, there can be additional solutions other than the solution studied above, 
namely, there are different regions where either V or / dominates. These new solutions 
can be much different from the minimal inflation models. In these models, a blue 
power spectrum is common and a negative running of the spectral index is possible. 
Large gravitational waves can be produced in these models, while in a large parameter 
region r lies below the experimental upper bound. Models of this kind with the 
quadratic and quartic potential are drawn in Figj6] and FigJTl 

FiglE] presents the /c-dependence of n s — 1 and a s for this solution. For the model 
we are considering, the running is not large enough for the blue spectrum to run to 
a red spectrum. 

By comparing with the n 3 = minimal models and the corresponding figures for 
the models with n 3 = 1, we find that a large running is usually produced in n 3 = 2 
models. While in the n 3 = 1 and n 3 = case, it is difficult to produce such a large 
running except in the small field quadratic model. 

It can be shown that the class of models / = M 2 R + AoM~ 4 yj 4 i? 2 are quite similar 
to models / = M 2 R + \ M~ 2 (p 2 R 2 , the only important difference is that in order to 
get the same n s — 1 and the running, Aq required in the former models are commonly 
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Figure 7: tit and r as functions of Ao in the models with / = M 2 R + XoM~ 2 ip 2 R 2 for the 
new solution in the nonminimal model. Tensor mode perturbation is large in this solution, 
which may be approved or disapproved in the near future. 




Figure 8: n s — 1 and a s as functions of ln(k/k c ) in the models with / = M 2 R+ XqM~ 2 ip 2 R 2 
and the quadratic potential V(ip) = ^m 2 ip 2 when / dominates over V. Ao is taken -10000 
in both figures. In order that r is not too large, the running can not be large enough to 
change the blue spectrum to a red one. 
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Figure 9: n s — 1 and a s as functions of Ao in the the models with f = M^R + 
Ao exp(— Xiip)R 2 . A blue spectrum is produced in the quadratic potential, but the run- 
ning is positive. A number —1.00E + 014 in the figure denotes —1.00 x 10 14 . 

3 or 4 magnitudes smaller than that in the latter models. 

6.3 Models with / = M*R + A exp{- \i(p / M p ) R 2 

In this class of models, to make any differences from the minimal model, a very large 
Ao is required. A blue power spectrum can be produced with a large but positive 
running. n s — 1, a s , and n?, r are shown in FigJHland FigfTTJl respectively 

7 Conclusion 

We obtained the scalar power spectrum V-ji, the spectral index n s — 1, the running 
of the spectral index a s , the tensor mode spectral index ny, and the ratio of the 
amplitude of tensor fluctuations to scalar fluctuations r in this paper, they are given 
in equations fpf3l) . (jHJ), (H61) . (jUJ) and fH9l) respectively. We constructed a nearly- 
conserved quantity which can be used in models different from the single inflaton 
model with minimal coupling. It has a wide range of applicability. Also, a new 
consistency relation is given, which is subject to test in future experiments. 

In the study of some concrete models, we see that there can be novel new features, 
such as a blue spectrum and a large running. We can also arrange F properly in order 
to have a brief fast-rolling before the slow-rolling period. It has been shown that the 
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Figure 10: n<r and r as functions of Aq in the models with / = M^R + Aq exp(— \\ip)R 2 . 



problem of a large running of the spectral index is eased in our models in contrast to 
the minimal model or f((p)R coupling models. 

The following issues deserve further investigation: The case when a nonminimal 
coupling becomes so strong that the quantum gravity effect destroys the usual canon- 
ical quantization picture at a length scale much smaller than the Hubble radius; the 
region where F crosses zero, producing a sudden fast-roll stage during the period of 
slow-roll inflation. A more careful study of these issues may lead to very interesting 
physics. 
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A A General Nearly Conserved Quantity 

From a general differential equation with two small parameters e\ and 62, 




(52) 
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The following quantity is nearly conserved when q <C aH, 

K = c[^i + (1 + e x + JL - e 2 )<j>)ef^ Hdt (53) 
where c is a constant. 71/ (H1Z) is of order 0(e 2 ) plus a term proportional to q 2 /a 2 . 



H v „ / ff 



,2 



= c[6 2 ( ei + - e 2 )if + ei + I I - eafoe/*** - (54) 

This construction is unique up to the constant c, and can be used in various cases 
where ( 152]) is satisfied. 

In the minimal inflation model, the equation (152]) takes the form 

} + (ff - | )0 + (2if - ^)0 + ^0 = (55) 
ii H a 

And 7?- is the comoving curvature perturbation, which happens to be conserved to all 
orders in e. 

In our nonminimal model, Inserting the Taylor expansion of Qk into (1251) . one can 
find 

<t>"k + (aH) L (1 2 "° Go) + M + 37 j & + A 

(56) 

where q is a rescale of k as in (1291) . So 7?. can be written as 

H ■ f H 2 FH G H \ 
K ^ #i+ ^ + 2 __ ? + __ f _ 1 j^ (57) 

It is shown in Appendix B that calculations using the comoving curvature per- 
turbation can also get correct results in our nonminimal model. But it can be shown 
that this fact is not generic in other generalizations of inflation models. 

The only requirement for 7Z to be nearly conserved is a diagonized equation of 
motion. So it can also be used in many other cases such as inflation models with 
holographic dark energy components [17J where entropy perturbation is produced 
and the comoving curvature perturbation is not a conserved quantity. 
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B Entropy-Perturbation-Like Sources 



In the minimal inflation model, during inflation and the post-inflationary period be- 
fore the comoving wave length re-enter the horizon, the pressure can be expressed 
as p = p(p, S), where p is the energy density and S is the entropy density up to a 
normalization. If there is no entropy perturbation produced, just as in the single- 
field minimal inflation model, there is only one degree of freedom in the scalar type 
perturbations, so p only depends on p. Then one can find 

" = T W 

p op 



This relation has two consequences. First, when (loo)) is satisfied, the equations 
f|T9|) and (l20l) can be diagonized into a one-variable differential equation exactly, as in 
the f(R) or f((p)R gravity models, studied in [T3]. Second, the comoving curvature 
perturbation is conserved in a model without entropy perturbation. 

But in general, when entropy perturbation SS is generated, the entropy perturba- 
tion can be defined as 

5p = ^5p + T5S (59) 
P 

where T is the temperature. It can be shown that 

~ H T k 2 
TZ = — ; — rj^SS + terms proportional to— r- (60) 

where TZ is the comoving curvature perturbation. 
During inflation, TZ can be expressed as 



a. 



H 

-R = 1p- -r- U + Hcf)) (61) 



H 

In nonminimal models, for general f((p,R), it is different from TZ introduced in 
Appendix A. From the uniqueness of TZ, we conclude that SS is nonzero. The existence 
of 5S can also be shown by direct calculations. 

It is not surprising that the entropy perturbation is produced in our nonminimal 
models, because when we consider nonminimal models other than f(R) or f(tp)R, one 
extra low energy effective degree of freedom opens up, and the entropy perturbation 
can be produced. 
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Although the entropy perturbation in our models exists and affects the evolution of 
scalar type perturbations, it can not live so long to have direct observable effects. To 
see this, we consider the period begin from a few e-folds outside the horizon, where the 
long wave length expansion can be applied to the Bessel function solution. During 
this period, ip and can be considered as slow-roll quantities, and the difference 
between time derivatives of 1Z and 1Z can be neglected. So 



K 



0(e 2 



(62) 



which means the entropy perturbation can be neglected a few e-folds after horizon 
crossing, and do not destroy the observational bounds. No extra decay mechanism 
for SS is needed. 
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